Change of ^Be decay rate in exohedral and endohedral Ceo fullerene compounds and 

its implications 



X 



A. Ray°, P. Das*^, S. K. Saha^ S. K. Das^ J. J. Das^ N. 
Madhavan^ S. Nath^ P. Sugathan^ P. V. M. Rao^-'', A. Jhingan^ 
'^Variable Energy Cyclotron Centre, 1/AF, Bidhannagar, Kolkata- 700064, India 
^ Radiochemistry Division, Variable Energy Cyclotron Centre, 1/AF, Bidhannagar, Kolkata- 700064, India 

'^Nuclear Science Centre, New Delhi, India and 
''Department of Nuclear Physics, Andhra University, Visakhapatnam - 530003, India 

(Dated: September 15, 2005) 

The half-life of ^Be implanted in a Ceo fullerene pellet and gold foil has been measured to be about 
the same within ~0.2%. Using a radiochemical technique, we also measured that the probability of 
. . . formation of endohedral ^BeQCeo by nuclear implantation technique was (5.6±0.45)%. It is known 

' from earlier works that the half-life of endohedral ^BeOCeo is ~1.2% shorter than that of ^Be 

' implanted in gold. An analysis of these results using linear mufRn-tin orbital method calculations 

, indicates that most of the implanted ^Be ions in fullerene Ceo stay at a distance of «5.3 A from the 

(N . 

centers of nearest Coo molecules forming exohedral compounds and those who enter the fullerene 
^ ^' cages go to the centers of the cages forming endohedral ^BeQCso compounds. 

(D ' 

. PACS numbers: 21.10.Tg,23.40.Hc,61.48.-|-c,23.40.-s,36.40.-c 

in 

It is known from earlier works0, |^ ^ IM IS that the rate of orbital electron capture of ''Be is susceptible to 
the surrounding environment and depends on both the lattice structure and electron affinity of the host atoms. 
Calculations 0, Isl, 'i^l also showed that the decay rate of ^Be should depend on its position in the host lattice. So the 
decay rate of ^Be in an atomic cluster such as fullerene Cgo or large biomolecule should also depend on its position 
with respect to the molecule and could be used as a tool to learn about the position of the implanted radioactive ion 
in the atomic cluster. In the future, this kind of study might also tell us about any abnormal change of DNA molecule 
, in a living cell. 

ly-j ■ It is already well-known[3,IM0| that different types of metal atoms(Be, Kr, Xe etc.) can be inserted into Cm fullerene 
' cage forming endohedral compound by nuclear implantation technique. Many theoretical studiesflfl llll Il2 . Q| 

■ have been done regarding the charge and equilibrium positions of the implanted ion in both endohedral and exohedral 
\ fullerene complexes, but there is no corresponding experimental measurement. So it is important to address these 

■ questions experimentally. By comparing the measured half-lives of exohedral ^Be-Cgo and endohedral ^BeQCgo with 
O that of implanted ^Be in another well understood material such as gold, we can learn about the charge and equilibrium 
5 positions of ^Be in endohedral and exohedral fullerene Cgo complexes. 

In this work, we have measured the difference of half-lives of ^Be implanted in a gold (Au) foil and fullerene (Cgo) 
pellet. A 25/im thick gold(Au) foil and a 500/im thick fullerene(C6o) pellet were bombarded by a 18 MeV ^Be beam 
from Nuclear Science Center, New Delhi, India. The ranges of 18 MeV '^Be in Au and Ceo are 12 /im and 40/im 
respectively. The intensity of the ^Be beam was about 40,000 particles/sec and the duration of the irradiation 24 hours 
Ci I for each sample. In order to obtain a 18 MeV ^Be beam, a liquid nitrogen cooled hydrogen gas cell at one atmospheric 
pressure was bombarded with a 15 pnA 21 MeV ^Li beam obtained from the pelletron machine of Nuclear Science 
Center, New Delhi, India. The ^Be nuclei produced by the reaction ^ H{J Li,'^ BeYn at 0° were separated from the 
primary ^Li beam by using a recoil mass spectrometer called Heavy Ion Reaction Analyzer (HIRA) 15] operated in a 
suitable ion optical modejlg. Primary ''Li beam particles were rejected by a slit system installed at the intermediate 
focal plane of the spectrometer. The primary beam rejection factor was about 10^^ and the purity of the separated 
^Be beam was about 92%. The advantages of rejecting the primary ^Li beam and using only a high purity ^Be beam 
for implantation are to minimize the radiation damage of the sample and avoid production of unwanted radioactivity. 

The ^Be implanted samples were brought to Variable Energy Cyclotron Center, Kolkata, India for off-line counting. 
Following electron capture, a ^Be nucleus has a 10.4% probabihty[l3 of populating the first excited state of ^Li 
which decays subsequently to its ground state emitting a 478-keV 7-ray photon. The half-life of ^Be was determined 
by monitoring the intensity of this 478 keV 7-ray line with time. Two HPGe detectors (detector-1 and detector-2) 
having efficiency «30% were used to count the samples (^Be implanted Au foil and ^Be implanted fullerene pellet). In 
addition, a standard ^'^'^Ba source was also placed in front of each HPGe detector. The detectors were well shielded by 
lead bricks to avoid any cross-talk between them and also to reduce the background level. A typical 7-ray spectrum 
(for ^Be in Cgo fullerene) as recorded by detector-1 is shown in Fig.^ Apart from 478 keV 7-ray line coming from 
the decay of ''Be and other 7-ray lines from ^•^'^Ba, we also see standard background 7-rays such as 239 keV, 511 keV, 
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583 keV, 609 keV, 727 keV, 911 keV lines, but there was no other contaminant. The 7-ray spectrum of ^Be in Au is 
very similar to Fig. 1. 

The count rate of 478 keV 7-ray was about 1 count per sec at the beginning of the run. Both the HPGe detectors 
were started at the same time, data was accumulated for 24 hours, stored in a computer and the spectra were cleared 
and the counting restarted. After counting for 7 days, the positions of the samples were interchanged and counted 
again. This was done to take care of any systematic error. The counting was continued for about 6 months. 

From each day's spectra, we determined the counts under 478 keV (N^(478)) and 356 keV (N^(356))photo-peaks 
coming from ''Be and ^^'^Ba respectively. The ratio (478) /N^ (356) should be independent of computer dead time 
and systematic errors and decay exponentially with time. In the case of data sets for ^Be in fullerene Ceo and ^Be in 
Au samples taken with detector-1, the reduced chisquare values of the exponential fits are 1.5 and 1.6 respectively. We 
also have similar data taken by detector-2 for ^Be in fullerene Ceo and Au and the reduced chisquare values of these 
exponential fits are 2.7 and 2.9 respectively. The 7-ray line shape of detector-2 showed slight tailings on both sides and 
this is probably responsible for obtaining comparatively poor exponential fit from the data set recorded by detector-2. 
So we compared data sets taken by the same detector and analysed the same way with similar quality of fits (reduced 
chisquare values) to cancel out the effect of any systematic error. We obtain the percentage difference of the two 
decay rates i.e aU^)'^""^ =(0.060 ± 0.405)% and (0.087 ±0.264)% for detector-1 and detector-2 respectively, where 
only statistical errors have been considered. Taking weighted average, we finally obtain the percentage difference in 
decay rates of ^Be in gold and fullerene to be ^^^^^^^^^^^^=(0.079 ± 0.221)%. So we find that the half-life of ^Be in 
gold and fullerene is about the same within wO.2%. 

In order to determine the half-hfe of ""Be in Au, we used known half-hfe Tt^ of ^^^Ba =(3836±15)days and only the 
data set taken by detector-1. The absolute half- life of ^Be in Au as obtained from the data set (reduced chisquare 
value = 1.6 for exponential fit) taken by detector-1 is = (53.60±0.19)days, where only statistical errors have been 
considered. The half-life of ^Be in Au as obtained from the data set taken by detector-2 agrees very well with that 
obtained from detector-1. However we think it would not be appropriate to reduce error bar on the absolute value 
of half-life by combining results from two data sets which give significantly different reduced chisquare values when 
fitted with exponential functions. Norman et al.[i( obtained half-life of ''Be in Au = (53.311±0.041)days and the 
reduced chisquare value for exponential fit of their data set was = 1.04. However in their experiment, a primary ^Li 
beam was incident on a Au foil and so some radiation damage of the Au lattice was expected. We also performed a 
similar experiment by bombarding a kapton foil with a 48 MeV ^Li beam from BARC-TIFR pelletron machine and 
implanting recoiled ^Be ions emitted at 0° in a Au foil placed behind the target. The implanted sample along with a 
-'^^'^Ba source were counted using a HPGe detector for about two months and the half-life of ^Be was determined by 
monitoring the ratio of counts under 478 keV and 356 keV 7-ray photo-peaks with time. The reduced chisquare of 
the exponential fit was — 1.45 and taking the half-life|0| of -'^^'^Ba =(3836±15)days, we obtained the half-life of ''Be 
in Au = (53.328±0.082)days, in good agreement with Norman et al.'s j3|result. So we think the radiation damage 
of Au lattice caused by a high flux incident primary heavy ion beam might be resDonsible|l8jfor slightly lowering 
(0.54±0.36)% the half-life of ^Be in Au. 

In a subsequent experiment, we performed a radiochemical separation of endohedral ^BeQCeo compound. ^Be 
implanted fullerene Ceo catcher was dissolved in 5 ml 1,2,4 trichlorobenzene. Equal volume of 6N hydrochloric acid 
containing carriers was added to the solution and the mixture was thoroughly shaken for 5 minutes in a separating 
funnel. The organic and the aqueous phases were separated and the activities of both the phases were determined 
accurately by 7-spectroscopy using a 20% HPGe detector. The aqueous fraction was again mixed with equal volume of 
organic solvent and a second extraction was carried out. The second organic fraction contained <4% of radioactivity 
compared to the first extraction. The organic fraction was also filtered through a millipore filter (pore size = 0.45 mm) 
to remove any insoluble material. No activity could be detected in the filter paper fraction. It is known'7|from Ohtsuki 
et al.'s high pressure liquid chromatography (HPLC) work that the ^Be radioactivities observed in the organic phase 
come from molecules having similar mobility as Ceo molecules implying ''Be is somehow attached to fullerene Ceo 
molecule. So the possibilities are ^Be activities in organic phase are coming from either endohedral ^Be@Ceo complex 
or heterofuUerene complex or exohedral ^Be-Ceo complex. Since only elements of group IVb and Vb of periodic table 
are knownQ to form heterofuUerene compounds, so beryllium is not expected to form a heterofuUerene compound 
with Ceo ■ ''Be ions stopped in the interstitial positions forming exohedral complex with Ceo fullerene are not expected 
to form strong covalent bond with carbon atoms of Ceo and should be readily dissolved in the hydrochloric acid and 
remain in the aqueous fraction. Therefore, as concluded earlier0 also, the ^Be activities observed in the organic phase 
should be associated only with the formation of endohedral ''Be@Ceo fullerene complex. So the yield of radioactive 
endofuUerene can be accurately determined from the ratio of activities present in the organic and aqueous fractions. 
We found that the yield of endohedral ^BeQCeo i-e. the probability of insertion of ''Be in Ceo by nuclear implantation 
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technique is (5.6±0.45)%, where the estimated uncertainty includes both the statistical and systematic errors. 

The radiochemically separated organic fraction containing endohedral ^Be@C6o complex was dried in a plastic 
crucible. Two HPGe detectors having efficiency «30% were used to count endohedral ''Be@C6o and ""Be in Au 
samples along with a ^-^^Ba source for 3 months following identical procedure as described earlier. However the 
amount of extracted endohedral ''BeQCeo complex was only 5.6% and so the statistical uncertainty on one day count 
of 478 keV 7-ray from endohedral complex was rather large (initially about 4.5% and 14% at the end of the run) for 
each detector. We finally obtained from our measurement that the half-life of endohedral ^BeQCeo complex is shorter 
than that of ^Be in Au by (3.3±2.3)%. Recently Ohtsuki et al.^^ measured half-life of endohedral ^BeQCgo complex 
with high accuracy and found that the half-life of endohedral ^BeHCeo complex is shorter than that of ^Be in Au by 
(1.20±0.12)%, in agreement with our measured half-life difference(3.3±2.3)% within one standard deviation. 

Let us first try to get a qualitative understanding of our results in terms of electron affinity. We have found 
experimentally that an implanted ^Be ion has a very low probability (5.6±0.45)% of entering the fullerene cage 
and they mostly stay in interstitial space forming exohedral ^Be-Cgo complex. The electron affinities of a fullerene 
molecule(C6o) and gold atom are 2.6 eV j2(| and 2.3 eV respectively. Both gold and fullerene have face-centered 
cubic (FCC) lattice structure, but the lattice parameter of fullerene Ceo lattice is much larger(14.17 A) than that(4.08 
A) of Au lattice. So if an implanted ^Be ion occupies the same geometrical position in both Cgo and Au lattices, then 
its distance from the center of the nearest Au atom and Ceo molecule would be very different. When ^Be would occupy 
the octahedral site of Au lattice, then its distance from the nearest Au atom would be Ri2A, whereas in the case of 
Ceo lattice, the corresponding distance would be «7A and the distance from the nearest carbon atom of Ceo molecule 
would be «3.7A. Hence ''Be should retain a significantly higher fraction of its 2s electrons in Ceo lattice and so the 
decay rate of ^Be should be significantly faster in Ceo lattice than in Au lattice. So the observation of the same(within 
0.2%) half-life of ^Be in Ceo and Au lattices should imply different geometrical positions for the implanted ''Be ions 
in the two lattices. It should also imply that the bond lengths of ^Be-Cgo (i.e. the distance between ^Be and the 
nearest carbon atom of Ceo) and ''Be-Au are about equal. 

In order to understand the results quantitatively, we have done tight binding linear muffin-tin orbital(TB-LMTO) 
method calculations IH 1^ IHi to determine the average number of 2s electrons of ^Be for both endohedral and 
exohedral fullerene complexes. A fullerene molecule consists of 60 carbon atoms placed on a sphere of radius 3.54 A. 
These fullerene Cgo molecules are arranged in a face-centered cubic structure with 14.17 A lattice constant. For both 
endohedral and exohedral complex, ''Be should go to equilibrium positions where the total energy of the system has 
local minima and the calculated decay rate of ^Be at those positions should agree with the experimental results. 

However we found that TB-LMTO code 22] was not very suitable for the determination of equilibrium positions in 
the fullerene Ceo lattice, because its spheridization of potential requires filling up of the available empty space with 
close packed empty spheres and since such division is not unique, it introduces an uncertainty in the determination 
of the total energy of the system. In the case of fullerene Ceo lattice, this uncertainty could be of the order of 0.4% in 
some cases and this is unacceptably large because the differences in binding energies of ^Be at different sites are very 
small (of the order of 10 eV). The uncertainty in the determination of the charge density and number of 2s electrons 
of ^Be could be up to ten percent in some cases, but this is acceptable for the purpose of estimating the change 
in half- life of ''Be and comparing with the experimental results. So the equilibrium position of ^Be in endohedral 
^Be@Ceo was taken from a previous density functional calculation 13J and TB-LMTO code was used to determine 
the average number of 2s electrons of ''Be at that position. 

Lu et al.|l3l| performed density functional calculations for endohedral Be@Ceo and found that the equilibrium 
position of beryllium ion should be exactly at the center of Ceo cage. So we performed TB-LMTO calculation 
placing a ^Be atom atom exactly at the center of fullerene Ceo cage. The calculation was performed assuming such 
endofuUerenes will form a face-centered cubic lattice with lattice parameter 14.17 A. Empty spheres were placed in 
the interstitial space and one inside the cage for close packing. TB-LMTO calculation shows that the radius of the 
muffin-tin sphere of the ^Be atom at the center of Ceo cage is «3 A and it almost fills up the entire space inside the 
Ceo cage. Let "^totai be the complete electronic wave function and ^'_Be2s be beryllium 2s state wave function. Then 
the square of the overlap of '^totai with '^Be2s i-c | < ^tota;|^Be2s > |^ represents the average number of 2s electrons 
in beryllium ion. Our TB-LMTO calculation shows that the average number of 2s valence electrons of ^Be atom 
placed at the center of a Ceo cage is — 1.07. 

We also calculated the average number of 2s electrons of ^Be in Au. As shown earlier(g],in the case of implantation 
of ^Be in Au lattice, ^Be should go to octahedral and tetrahedral sites of Au lattice. TB-LMTO calculations were 
performed placing a ^Be atom at octahedral and tetrahedral positions of a face- centered cubic Au lattice having 
lattice parameter 4.08 A. The average number of 2s electrons | < ^'tota;|^Be2s > |^ was found to be 0.54 and 0.36 for 
^Be occupying octahedral and tetrahedral positions of gold lattice respectively. We expect that as a result of random 
implantations, the number of ^Be atoms occupying tetrahedral sites would be twice that of occupying octahedral 
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sites, because the number of tetrahedral sites are twice that of octahedral sites in a face-centered cubic lattice. Hence 
taking weighted average, we found ^Be would have 0.42 electrons in its 2s orbital when implanted in Au. So the decay 
rate of ^Be in endohedral ^BeCSCgo should be faster than that of '^Be in Au. According to Hartree and Hartree's 
calculation [2^, the overlap of valence 2s electrons to the total electronic overlap at beryllium nucleus is only 3.2% 
for a neutral beryllium atom. As shown in ref[6], there is a linear relationship between the decay rate of ^Be and 
its number of valence 2s electrons. Using Hartree and Hartree's result^^ along with the linear relationship found 
earlierj^, we obtain that the decay rate of endohedral '"BeQCgo should be faster by ~ 1.1% compared to that of ''Be 
in Au lattice. This result agrees well with our measurement and Ohtsuki et al.'s earlier measurement [l that the 
decay rate was faster by (1.20±0.12)%. 

We have found that the measured probability of formation of endohedral ^BeQCgo complex by nuclear implantation 
technique is only 5.6% and the remaining 94.4% of the time, ^Be goes to the interstitial space forming exohedral ^Be- 
Cgo complex. In order to calculate the decay rate of ^Be when it goes to interstitial space of Cgo lattice, we have to 
first determine its equilibrium positions in the fuUerene Cgo lattice. 

Although we have not come across any calculation for determining equilibrium positions of ^Be in a fuUerene Ceo 
lattice, there are a large number of calculations for equilibrium geometries of isolated exohedral fuUerene 
Ceo complex. From those calculations, it is known that the typical bond length between a carbon and the other ion 
formin g ex ohedral fuUerene Ceo complex is w 2 A. The distance of the other ion from the center of the Ceo cage is 
usuallv|lol Ull around 5.5 A and the equilibrium position of the other ion is most likely to be along a line 

bisecting normally a C-C bond of Ceo molecule. Using this information, we carried out TB-LMTO calculations for 
beryllium in Ceo lattice and found that there were total energy minima postions when beryllium was on any face of the 
face-centered cubic Ceo lattice between a corner Ceo molecule and a face-centered Ceo molecule along a line bisecting 
normally C-C bond of a Ceo molecule. At equilibrium position, the distance of beryllium ion from the center of Ceo 
molecule is about 5.3 A and C-Be bond length comes out around 2 A. We find from TB-LMTO calculation that at 
this position, the average number of valence 2s electrons of beryllium ion is = 0.44. So using Hartree and Hartree's 
result ^23*1 along with the experimental observation that only 5.6% of the implanted ^Be ions form endohedral complex, 
the decay rate of the implanted ^Be in fuUerene Ceo should be faster by w 0.1% compared to that of ^Be in Au. This 
result agrees with our experimental observation that the decay rates of ^Be in fuUerene Ceo and Au are equal within 
0.2%. 

If we assume that ^Be ions go to the octahedral or tetrahedral sites of face-centered cubic Ceo lattice, then according 
to TB-LMTO calculations, the average number of 2s electrons of ^Be should be 1.12 and 1.51 respectively. So using 
the method of rcf[6], the decay rate of "^Be should be faster than that of ^Be in Au by 1.16% and 1.8% respectively, 
in complete disagreement with our experimental observation. Moreover if ^Be is at octahedral site, then its distance 
from the center of the nearest Ceo molecule would be about 7 A and so the length of shortest C-Be bond would be 
about 3.7 A. Similarly when ^Be is at tetrahedral site of Ceo lattice, then its distance from the center of the nearest 
Ceo molecule is 6.1 A and the length of shortest C-Be bond is about 3 A. At this position, the line joining the center 
of the Ceo molecule and ^Be ion does not bisect C-C bond of Ceo molecule, but passes through the center of the 
hexagonal opening of Ceo molecule. It is also expectedjjjjj that for most stable exohedral configuration, a beryllium 
atom should go to a site of high electron density. Our TB-LMTO calculations show that the electron density at 
tetrahedral site of Ceo lattice should be about 50% more than that at octahedral site. However the electron density 
on a face of the lattice midway between two Ceo molecules is about five times the electron density at the octahedral 
site. So from these qualitative arguments, we expect that the octahedral and tetrahedral sites should not be the 
equilibrium positions of ^Be in fuUerene Ceo lattice. However our TB-LMTO calculations for the total energy of the 
system do not support this expectation, perhaps due to the model limitations mentioned above. 

In the case of endohedral ^Be@Ceo complex formation, although C-Be bond length would be 3.54 A since ^Be is 
expected to be at the center of fuUerene Ceo cage, however at that position, ^Be is equidistant from all 60 carbon ions 
and is held there by sixty C-Be bonds. So even if each bond is weak, there are sixty such bonds holding the beryllium 
ion. 

In summary, we have found experimentally that the decay rate of ^Be in fuUerene Ceo is equal to that of ^Be in Au 
within 0.2%. A radiochemical analysis of the irradiated sample shows that only 5.6% of the implanted ^Be ions form 
endohedral ^Be@Ceo complexes and the remaining 94.4% of the ^Be ions form exohedral ^Be-Ceo complexes. The 
measured decay rate of endohedral ^Be@Ceo complex is in agreement with recent accurate measurement of Ohtsuki 
et al.^^. Our TB-LMTO analysis shows that for endohedral ^Be@Ceo complex, ^Be ion should go to the center of 
Ceo cage. For exohedral ^Be-Ceo complex, we expect that ^Be should go to one of the faces of face-centered cubic 
Ceo lattice forming exohedral complex with C-Be bond length «2 A. 

We acknowledge useful discussion with R. Vandenbosch (University of Washington, USA) ,Ove Jepsen (Max Planck 
Institute, Germany) and J. Lu (Peking University, Beijing, China), G. P. Das (lACS, Kolkata, India) and L. M. 
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FIG. 1: 7-ray spectrum from the decay of '^Be implanted in Ceo fuUerene 

Ramaniah (Bhabha Atomic Research Centre, Mumbai, India) and A. Mookerjee (S. N. Boss Centre for basic sciences, 
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